Anomalous flux avalanche morphology in a a-MoGe superconducting fllm with a 
square antidot lattice - experiment and simulation 
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We have employed magneto-optical imaging to visualize the occurrence of flux avalanches in a 
superconducting film of a-MoGe. The specimen was decorated with square antidots arranged in a 
square lattice. We observed avalanches with the anomalous habit of forming trees where the trunk is 
perpendicular to the main axis of the square lattice, whereas the branches form angles of 45 degrees. 
The overall features of the avalanches, and in particular the 45 degree direction of the branches, 
were confirmed by numerical simulations. 

PACS numbers: 



Flux avalanches triggered by thermomagnetic instabil- 
ities have been reporte d in a variety of superconduct- 
ing films of Nbi, MgB.^, NbgSiP, NblSf^, YNi2B2C^ 
and YBCC0. Such events, in which flux bursts suddenly 
invade the sample, have also been studied in supercon- 
ducting films decorated with a lattice of antidots, which 
guide flux motion, as r evea led by magneto-optical imag- 
ing (MOI) experimenttPH!. 

The sample investigated in the present work is an 
amourphous film of MoGe (a-MoGe) approximately rect- 
angular in shape, with thickness t = 25 nm and lateral 
dimensions roughly 2.0 mm x 2.6 mm. A square lat- 
tice of square antidots (ADs) was fabricated by electron 
beam lithography. The a-MoGe was deposited via pulsed 
laser deposition from a target of Mo78Ge22 with purity 
99.96%, on top of a Si02 insulating substrate. The pat- 
terned sample has square antidots with sides 0.4 /xm. The 
period of the pattern is w = 1.5 fim, which corresponds 
to a commensurability field i/i = ^ = 9.2 Oe at which 
the densities of vortices and antidots match each other^^. 
Using the expression for the upper critical field and for 
the dirty limilfi^, the zero temperature superconducting 
coherence length, ^gl(O) = 5 nm was determined, and 
the penetration depth, Xcl{0) = 517 nm was obtained 
for a similar plain film. Figure lb shows the orientation 
of the antidot lattice relative to the edges of the sample. 

Characterization measurements of the ac susceptibil- 
ity and dc magnetization were carried out in commercial 
Quantum Design equipments (MPMS and PPMS) show- 
ing that the superconducting transition temperature at 
zero field is Tc = 6.74 K. The MOI technique employed 
relies on the Faraday effect a garnet indicator fllm 
placed on top of the superconducting specimen. The in- 
dicator used in the present work was a Bi-substituted 
yttrium iron garnet film (BiiYIG) with in-plane magne- 
tization. 



Figure la shows a magneto optical image taken at T = 
4.5 K and an applied DC field of H = 1 Oe. Noticeably, 
the main trunk of the tree-like avalanches are predomi- 
nantly perpendicular to one of the sample edges whereas, 
surprisingly, the branches form angles of 45 degrees rel- 
ative to the main axes of the antidot lattice. Avalanches 
directed at 45 degrees have previously been observed only 
in a Nb film and there only from edges strongly inclined 
relative to the main axis of the antidot lattice^. 

Figure 2 shows the results of three different experi- 
ments, all of them carried out at T = 3 K with H = 1.4 
Oe. The images are colored red, green and blue, so that 
overlapping penetrations in all three experiments appear 
as shades of gray, whereas non-repetitive flux patterns 
combine in color. Evidently, there is essentially no over- 
lap in the dendritic flux patterns formed in the three 
experiments. The background flux penetration, however, 
is fully reproducible. 

The origin of dendritic avalanches in superconducting 
fllms is a thermomagnetic instability mechanism between 
the Joule heating created by vortex motion and the re- 
duction of the critical current density as temperature in- 
creases. The instability is also a consequence of the non- 
linear material characteristics of type II superconductors, 
which is conventionally approximated by a power law 



where E is electric held, J is sheet current, J = |J|, 
po is a resistivity constant, d is sample thickness, Jc is 
critical sheet current, and n is the creep exponent. The 
temperature dependencies are taken as 

Jc^ Jcoil-T/T,), n^m/T, (2) 

where Tc is the critical temperature. The electrodynam- 
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FIG. 1: (a) Magneto optical image of flux penetration at T 
= 4.5 K and H = 1 Oe, applied perpendicular to the film. 
The image brightness represents the magnitude of the local 
flux density, (b) Optical image of the square lattice of square 
antidots. 



ics must be supplemented by the heat diffusion equation 
cf ^ n\7^T-'^{T-To) + ^JE, (3) 

where c is specific heat, k is thermal conductivity, h is 
the coefficient for heat removal to the substrate, and Tq 
is the substrate temperature. 

Eq. ([3| must be solved together with Maxwell's equa- 
tions and the material law, Eq. ([T]). The description of 
the simulation method is found in Ref. [TH so here only 
the key ingredients are outlined. The main challenge is 
to invert the Biot-Savart law in an efficient way. This is 
done by including also the vacuum outside the sample in 
the simulation formalism. At the cost of including the 
extra space, one can use a real space/Fourier space hy- 
brid method with the very attractive performance scaling 
of 0{N log N), where N is the number of grid points. 

The area of the sample was discretized on a 512x512 
grid, with the geometry of strip, emulated by periodic 
boundary conductions. A total number of 43x43 square 
antidots where imposed, so that the antidot coverage was 
1 /4 of the total area. The initial state of the simulation 
was prepeared by ramping the external field with thermal 
feedback turned off, so that a critical state was formed 
from the edges. When later, the thermal feedback was 
turned on, and at the same time a local spot close to the 
edge was heated above Tc, a flux avalanche developed 
from the heated spot. 

Figure 3 shows a map of the local flux density obtained 
by the numerical simulation. The avalanche of the figure 
demonstrates a pronounced guidance effect, contrary to 
the case of plain films or films patterned with random 
disorder. Just as in the experiments, there is a main 





FIG. 2: Overlay of 3 repeated identical MOI experiments 
showing the flux penetration into the a-MoGe film, carried 
out at T = 3 K with H = 1.4 Oe. Each experiment is color 
coded in red, green and blue, thus the presence of these colors 
in the overlay shows that the penetration is non-reproducible. 



FIG. 3: Results of numerical calculations of the flux penetra- 
tion into a superconducting film with a square array of square 
antidots where a thermomagnetic feedback is taken into ac- 
count. Details of the numerial procedure will be published 
elsewhere. 
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trunk perpendicular to the edge and two main branches 
forming an angle of 45 degrees with both the edge and 
the directions of the antidot lattice. 

In summary, we have employed MOI to observe flux 
avalanches in superconducting films of a-MoGe decorated 
with a square lattice of square antidots. Avalanches have 
the form of trees where the main trunk is perpendicular 
to the sample edge, whereas its branches form an angle 
of 45 degrees with the main axes of the antidot lattice. 
The overall features of the avalanches, and in particular 
the 45 degree direction of the branches, were confirmed 



by numerical simulations. 
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